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We determine the adiabatic phase diagram of a resonantly-coupled system of Fermi atoms 
and Bose molecules confined in the harmonic trap by using the local density approximation. 
The adiabatic phase diagram shows the fermionic condensate fraction composed of condensed 
molecules and Cooper pair atoms. The key idea of our work is conservation of entropy through 
the adiabatic process, extending the study of Williams et al. [Williams et al, New J. Phys. 6, 
123 (2004)] for an ideal gas mixture to include the resonant interaction in a mean-field theory. 
We also calculate the molecular conversion efficiency as a function of initial temperature. 
Our work helps to understand recent experiments on the BCS-BEC crossover, in terms of 
the initial temperature measured before a sweep of the magnetic field. 

KEYWORDS: degenerate Fermi gas, Feshbach resonance, BCS-BEC crossover, adiabatic 
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1. Introduction 

Degenerate Fermi gases with controllable interaction have been realized in recent exper- 
iments by making use of a Feshbach resonance. Much of the current study of these systems 
focuses on the crossover between a Bardeen-Cooper-Schrieffer (BCS) superfluid of fermionic 
atoms and Bose-Einstein condensation (BEC) of diatomic molecules. 1-4 In these systems, the 
position of a molecular bound state can be tuned relative to the scattering continuum by vary- 
ing an applied magnetic field. By ramping the magnetic field across a Feshbach resonance, 
one can continuously transform a fermionic atomic gas to a bosonic molecular gas, both above 
and below the superfluid transition temperature. 

Several theoretical papers reported the phase diagrams of the resonantly-coupled Fermi- 
Bose mixture gas. 5-7 In order to compare these works with the experimental results by Regal 
et at 1 and Zwierlein et aZ., 3 one has to understand how the system traverses the phase diagram 
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as the resonance energy is varied. However, the experiments only provided a measure of the 
initial temperature, Tj, before the magnetic field sweep, thus complicating the attempts to 
relate the data to theoretical models. In the experiments, the sweep is slow so that atoms and 
molecules can move and collide sufficiently in the trap, 3 allowing for relaxation to thermal 
and chemical equilibrium during the sweep. Assuming the sweep to be adiabatic , Williams 
et al. s calculated the phase diagram for the ideal gas mixture of fermionic atoms and bosonic 
molecules as a function of the resonance energy and Tj. The key idea of this adiabatic process 
is the conservation of entropy during the sweep: the system follows a path of constant entropy 
through the phase diagram, and the entropy is set by the initial preparation of the sample, 
i.e. the initial temperature. Using the reverse logic, Chen et al. 9 gave the idea of thermometry, 
by which the real temperature of the system after an adiabatic sweep may be deducted from 
the initial temperature before the sweep. Based on this finite temperature formalism, they 
compared the experimental phase diagram with their theoretical boundary between the normal 
phase and the superfluid phase. 10 Similar ideas were proposed by Hu et al. 11 and Carr et al. 12 

In this paper, we extend the work by Williams et al. s for an ideal gas mixture of fermionic 
atoms and bosonic molecules to include the resonant interaction using mean-field theory. Our 
mean-field treatment allows for the description of the BCS superfluidity due to the interaction, 
which was neglected in ref. 8. In order to introduce the effect of the harmonic trap, we apply 
the local density approximation (LDA). 

Calculating the entropy as a function of the temperature and the detuning, we show the 
paths of constant entropy traversed in the conventional phase diagrams as the detuning is 
varied adiabatically. On the basis of the adiabatic path, we determine the adiabatic phase 
diagrams against the detuning and entropy, labeled by the initial temperature of the Fermi 
gas measured before a sweep of the magnetic field. We note that the adiabatic phase diagram 
of ref. 10 plots the superfluid fraction, whose phase boundary has a good agreement with the 
experiment in ref. 1, while the experimental phase diagrams 1,3 show the condensate fraction 
which is zero momentum molecules projected through the fast sweep, assuming that the 
fraction has an information of zero momentum molecules and zero center-of-mass momentum 
Cooper pairs. In this paper, we plot the phase diagrams for the condensate fraction that is 
composed of zero momentum molecules and zero momentum Cooper pairs. 

We are also interested in the production efficiency of molecules during the adiabatic sweep 
of the magnetic field. The conversion efficiency is found experimentally to be a monotonic func- 
tion of the initial peak phase space density, when the magnetic field is swept adiabatically. 13 
A recent paper by Williams et al. 1A discussed the mechanism of molecular formation, in the 
case where three-body recombination may be neglected, and derived an expression for the 
molecular conversion efficiency in terms of the initial peak phase-space density. Their result 
agrees with the experimental data 13 although they treat a non-interacting quantum gas. We 
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extend the calculation of the molecular conversion efficiency in ref. 14 to include the mean-field 
effect. 

Our study gives an intuitive understanding of recent experiments on the BCS-BEC 
crossover using the adiabatic sweep process. 

2. Equilibrium Theory 

Our model is the coupled boson-fermion Hamiltonian with a two-component Fermi 

gas; 6 ' 15 ' 16 



H 



= f dx[$, (r)# a (r)^(r) + ft (r)H m (r)4>(r)} 

+k J dr[ft(r)^(r)^(r) + ^(r)^(r)^(r)], (1) 

where repeated indices a imply a sum over pseudospin states a = { j, J,}. The field operators of 
atoms and molecules are respectively denoted as V> CT (r) and 4>(r), which obey Fermi and Bose 
commutation relations. The single particle Hamiltonians for the atoms and molecules are given 
by H a (r) = -(h 2 /2m)V 2 + V a (r) and H m (r) = -(h 2 /2M)V 2 + V m (r) + 5 , respectively, where 
m is the atomic mass and M = 2m is the molecular mass. We introduce the harmonic trap in 
which the each pseudo-spin Fermi atoms feel the same potential, given by V a {r) = mco 2 r 2 /2. 
On the other hand, we assume that molecules feel the following potential: V m (r) = Muj 2 r 2 /2. 
This assumption is valid for the current experiments in optical dipole traps. 17 Two fermionic 
atoms are coupled to a bosonic molecule through the resonant interaction with the coupling 
constant k. The detuning 6(B) can be tuned by adjusting an external magnetic field B, since 
the bare atoms and molecules have different magnetic moments, and thus experience a relative 
linear Zeeman shift. As usual we take the zero of energy to be the energy of two separated 
atoms at rest at each value of B. We neglect the non-resonant atom-atom interaction, since 
the effective interaction for two atoms is dominated by the resonant contribution —k 2 /5 near 
the resonance. 

In order to impose the constraint that the number of particles be conserved, we introduce 
the chemical potential and deal with the grand canonical Hamiltonian H = H — fiN, where 
the total number operator of particles is defined as N = J dr[J24>i(r)'4> (7 (r) + 2$t(r)</>(r)]. 
The factor of two in the last term reflects that each molecules are assembled from two atoms. 
We regard the system as in chemical equilibrium between atoms and molecules. 8, 18 

In the superfluid phase, we can separate the molecular field operator into two parts: 

<Kr) = (0(r))+0(r) 

= d>(r) + <Kr). (2) 

Here, 3>(r) represents the molecular condensate wavefunction, which corresponds to the order 
parameter, and 4>(r) represents the non-condensate molecules. In the mean-field theory, one 
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leaves out the Cooper pairs outside the condensate, and the Hamiltonian H reduces to the 
BCS-type Hamiltonian 15 

H = Jdr{*%r)[H m (r)-2»]*(r)} + Jdr{ft(r)[H m (r)- 2 M ]0(r)} (3) 
+ f dr{^(r)[H a (r) - /i]i(r)} + « [ dr[**(r)^ T (r)^(r) + ^J(r)^}(r)*(r)], 



where fluctuation terms like (fMj)^ have been neglected. 

This type of quadratic Hamiltonian can be diagonalized by the standard Bogoliubov trans- 
formation. In order to discuss the thermodynamics of the system, we evaluate the grand 
canonical potential Q,. The effect of a harmonic trap in the grand canonical potential Q is 
included by the local density approximation, 5 which adopts the effect of the local potential 
by replacing the chemical potential \i with \i(y) = \i — V a (r). 

The grand canonical potential including the effect of the harmonic trap potential is given 

by 



n = |dr([5-2Mr)]|$(r)| 2 + i|^ln{l- e -^W(r)] } 
+ /^{^ a) -Mr)-^(r)]-|ln[l + e-^)]}). 



(4) 
(5) 



Here, = h 2 k 2 /2M is the kinetic energy of a molecule, and = h 2 k 2 /2m is the kinetic 
energy of an atom. As usual (3 = 1/kT, where T is the temperature and k is Boltzmann's 
constant. At this step, we have already assumed that the atoms and molecules are in thermal 
equilibrium. In other words, the atom temperature and the molecule temperature are the 
same. 8 ' 18 The local quasiparticle excitation energy is E^r) = \J [e^ — fi(r)] 2 + |A(r)| 2 , where 
the local gap energy is | A(r)| = «|<&(r)|. Within our mean-field theory, the excitation spectrum 
of bosonic molecules has the finite excitation gap S — 2/j,. However, in the superfluid phase, the 
symmetry breaking requires a gapless spectrum with the Bogoliubov phonon mode. This can 
be incorporated into the theory by including the effect of fluctuations, 19 which are neglected 
in our mean-field theory. The gapless excitations are important at very low temperatures. We 
will discuss this point in § 3. 

The local gap equation is obtained from 6Q/5\&(r)\ =0 , which leads to 



k 2 f dk 1 



S - 2 M (r) = - J ^ Ek(r) tanh 



(6) 



On the other hand, the total number of particles N tot follows from the relation N tot = 
—dQ,/dfi, which gives 

iVtot = 2jdrMr)\ 2 + 2j^f B (et ) -,r) 
drdk 



(2tt) 2 



+ 2 — — JF{E^, r) 



E k (r) £ k (r) 



(7) 
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Here, the Bose and Fermi local distribution functions are respectively /s(e^;r) and 
f F {E k ; r), which are given by /^(e^; r) = l/{exp [^(e^ + 5- 2p{r))\ - 1} and f F {E k ; r) = 
l/[exp (f3Ek(r)) + 1] .In the right-hand-side of eq. (7), the first term represents the twice 
number of the condensed molecules 2N mc . The second term represents the twice number of 
the non-condensed molecules 2N m . The last term is the number of atoms N a . Thus, eq. (7) 
can be written as 

JVtot = 2N m + N a 

= 2N mc + 2N m + N a . (8) 

We also calculate the condensed pair number N c , which is composed of the number of 
condensed molecules N mc and Cooper pairs N p : 

N c = N mc + N p 

= J drMr)\ 2 + J dndraK^Cn)^^))! 2 , (9) 

where the second term of eq. (9) describing the number of the condensed Cooper pair is 
the maximum eigenvalue of two-particle density matrix. 20 Within our mean-field theory, the 
explicit expression for the number of Cooper pairs N p is given by 21 ' 22 

drdk (\A(r)\ 2 



N p = /^{^tanh 2 ^!) 

J (2vr)3 W(r) [ 2 J J 



(10) 



We regard eqs. (6) and (7) as the simultaneous equations for the chemical potential p 
and the local gap |A(r)|, for a given temperature T and detuning 5. In this paper, we set the 
bare coupling strength to a = y/pn = 0.4ep) where p is the peak density and e-p is the Fermi 
energy of a pure gas of fermionic atoms, which corresponds to the narrow Feshbach resonance, 
which is given by p = {2me-p) 5 /(37r 2 7i 3 ) and £f = (3iVtot) 3 ^o- This weak coupling constant 
allows for a mean-field treatment. Although this value does not correspond to experimental 
situation of the broad Feshbach resonance, discussed in ref. 23 (see also discussion in § 3), 
we do not expect any qualitative differences in equilibrium phase diagrams between narrow 
and broad Feshbach resonances. To avoid the divergence of the integral in eq. (6), we impose 
the Gaussian cut-off; ex.p { — [e / (2ef)] 2 } , whose cutoff energy scale is chosen to be the Fermi 
energy. 5 

3. Adiabatic Phase Diagram 

In this section, we discuss the adiabatic phase diagram. We first briefly review the ex- 
perimental procedures of refs. 1 and 3. In these experiments, the ultracold two-component 
Fermi gas is initially prepared at a magnetic field far detuned from the resonance position, 
corresponding to a pure atomic gas. In this weakly interacting regime the temperature of the 
Fermi gas is measured using time-of-flight imaging. The magnetic field is then slowly lowered 
to the vicinity of the resonance to allow the atoms and molecules sufficient time to move and 
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collide in the trap. This process may satisfy the condition that the gas be able to collisionally 
relax to equilibrium. 8 We regard this process as an adiabatic quasistatic process, because the 
magnetic field is varied while keeping the thermal and chemical equilibrium from moment to 
moment without exchanging heat or particles with the environment. The system will thus 
follow a path of constant entropy in the phase diagram. 

The total entropy S^ot is the sum of the atom entropy S a and the molecule entropy S m . 
These are given by 

and 

bm ~ kB J (2vr)3 { e ^ +S - Mr)] _ 1 ln{1 j J- {U) 

Before showing contours of constant entropy, we plot the conventional phase diagrams of 
the condensed molecular fraction r\ mc = 2N mc /N tot and the condensed Cooper pair fraction 
r/p = 2N p /N tot separately against the temperature T and the detuning 5 in Figs. 1 (A) 
and 1 (B). From Fig. 1 (B), we find that the number of Cooper pairs at low temperatures is 
peaked above the resonance in the BCS regime, while its contribution to the total condensed 
fraction r\ c is very small. This is because we deal with the weak coupling (narrow Feshbach 
resonance). In the case of the strong coupling (broad Feshbach resonance), the Cooper pair 
contribution in the fermionic condensate fraction becomes dominant, and the single channel 
model describes effectively such systems in the vicinity of the resonance. 24 However, as shown 
in ref. 23, the conventional phase diagrams are almost the same for both cases, when plotted 
against (/cpa,,)" 1 instead of 5 (kp is the Fermi wavenumber and a s is the s-wave scattering 
length), ref. 23 also shows that the behaviors of the total number of bosons (including both 
Cooper pairs in the open channel and bare molecules of the closed channel) are almost the 
same for the two cases. We thus expect that the adiabatic phase diagram discussed below also 
describes the qualitative behavior in the broad resonance case. 

In Fig. 1 (C), we show contours of constant entropy in the S - T plane. We recognize those 
lines, as the paths that are traversed as the detuning S is swept adiabatically, starting from 
the right side of the resonance at a large value of the detuning 5. We see that the temperature 
increases as the detuning S is lowered adiabatically. This temperature increase on the path of 
constant entropy can be understood as follows: first, due to the conversion of pairs of atoms 
into molecules, the system loses degrees of freedom; second, the condensed molecules does not 
contribute to the entropy; third, the existence of Cooper pairs also suppresses the entropy. 
For these reasons, the system must then heat up in order to conserve the entropy. 

The relation between the initial temperature of the atomic gas and the final temperature 
of the molecular gas has been given for the ideal gas mixture of fermionic atoms and bosonic 
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Fig. 1. (Color online) (A) Conventional phase diagram of the condensed molecular fraction ry mc . 
(B) Conventional phase diagram of the condensed Cooper pair fraction r) p . Note the change in 
colorscale between panels (A) and (B). Each fraction in the conventional phase diagram is plotted 
against the temperature T and the detuning S, normalized by the Fermi energy. (C) Contours of 
constant entropy. The paths are plotted against the temperature T and the detuning 5, as the 
detuning is varied adiabatically. We set the coupling strength a — 0.4ep. 



molecules by Williams et al. s In the low-temperature limit T/Tp -C 1, one finds 

1/3 



r tt 2 - 


1/3 ^ 


[l2C(4)_ 





(13) 



while in the high-temperature limit, T/Tp S> 1, the relation becomes 



6 l/3 e 4/3 



(14) 



Tf " " \T F/ 

The above relations can be derived by connecting the initial entropy of a pure atomic gas 
and the final entropy of a pure molecular gas. Similar expressions have been derived by Carr 
et al. 25 in connection with the idea of cooling a Fermi gas by adiabatically increasing the 
magnetic field starting from a molecular gas. 

In Fig. 2, we plot the final temperature in a deeply BEC region and a moderately BEC 
region as a function of the initial temperature. The solid line is obtained by a numerical calcu- 
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Fig. 2. The final temperature versus the initial temperature where the final detuning is 5/ef = —100 
(solid line) and 5/e-p = —4 (dotted line). The dashed and dot-dashed lines represent the high and 
low temperature limit. We set the coupling strength a = 0.4ep. 

lation in a deeply BEC region where the detuning is 5 = — IOOcf- The dotted line is obtained 
by a numerical calculation in a moderately BEC region where the detuning is 5 = — 4ep- The 
dashed line represents the high-temperature limit for an ideal gas, while the dot-dashed line 
represents the low-temperature limit for an ideal gas. The low-temperature approximation 
agrees with the both of numerical solutions. The high temperature approximation does not 
agree with the numerical solution of a moderately BEC region, because atoms appear when 
the temperature is not sufficiently high and the concept of connecting the entropy of a pure 
atomic gas and that of a gas consisting entirely of molecules breaks down. On the other hand, 
the high temperature approximation agrees with a numerical solution of a deeply BEC region 
except for the high initial temperature regime where some atoms remain at the end of the 
sweep. These agreements are reasonable since our model assumes the narrow resonance case, 
leaving out the Cooper pairs outside the condensate. The system is equivalent to the ideal gas 
mixture for temperatures above T c . 

In Fig. 3, we plot the fermionic condensate fraction r\ c = 2N c /N tot against the detuning 5 
and the initial temperature of the atomic Fermi gas measured before a sweep of the magnetic 
field, as in the experimental phase diagrams of refs. 1 and 3. The initial temperature T; is 
found by equating the total entropy to the initial entropy Stot($, T) = Si(Ti). 

In Figs. 1 (A) and 1 (B), the transition temperature is T c ~ 0.53Tf in the BEC region 
where the detuning is 8 = —4ep- On the other hand, in Fig. 3, the transition temperature in 
terms of the initial temerature is Tj )C ~ 0.21Tp at 5 = — 4ep- This behavior is consistent with 
Williams et al. 8 for an ideal gas model, where the transition temperature is T c ~ 0.51Tp and 
the transition temperature in terms of the initial temerature is T^ c ~ 0.21Tp at 5 = — 4ep- 
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Vmc + f]p i s plotted against the initial temperature T| and the detuning S. 

The difference in the transition temperature between our model and the ideal gas model 
becomes smaller when we translate the conventional transition temperature to the transition 
temperature in terms of the initial temperature. This translated transition temperature Tj iC in 
the BEC limit is comparable with transition temperature observed in the experiments. 1 ' 3 Here 
we comment on the finite energy gap of bosonic excitations discussed in § 2. Because of this 
artificial gap, our mean-field results should be unreliable in the low temperature region where 
ksT < 5 — 2/i. However in our adiabatic phase diagram, for 5 < 0.5£F) the real temperature 
satisfies the inequality fc^T 3> 5—2fi as long as the initial temperature satisfies /ceTj > 0.05£f- 
The region where the real temperature satisfies ksT S> S — 2\i becomes wider as one goes into 
the BEC side of the resonance. On the other hand, for 5 > 0.5ep, the total number of molecules 
is so small at low temperatures, that the bosonic excitations make no significant contribution 
to the thermodynamic. Therefore, the finite energy gap in the bosonic excitation spectrum due 
to the mean-field treatment does not affect our adiabatic phase diagram significantly except 
for the low temperature region fceTj < 0.05£f in the crossover region 8 ~ 0. 

In Fig. 4, we plot the fermionic condensate, the condensed molecular, and the total molec- 
ular fractions with an initial temperature Tj/Tp = 0.08 as a function of the detuning. As the 
detuning is lowered, the fermionic condensate fraction increases. Subsequently it decreases 
after having a peak near the resonance. This behavior is consistent with experiments. 1 ' 3 
Williams et al. 8 finds a similar result for an ideal gas mixture, though they normalize the 
condensate fraction by the number of molecules. 

4. Molecular Conversion Efficiency 

Hodby et al. 13 have investigated the molecular conversion efficiency in bosonic 85 Rb and 
fermionic 40 K. They found that the molecular conversion efficiency does not reach 100% even 



9/13 



J. Phys. Soc. Jpn. Full Paper 




Fig. 4. (Color online) The fcrmionic condensate fraction rj c , the condensed molecular fraction rj mc and 
the molecular fraction rj m for the initial temperature Tj = 0.08Tp are plotted against the detuning 
6. The solid line represents the fermionic condensate fraction r/ c . The dotted line represents the 
condensed molecular fraction r\ mc . The dashed line represents the molecular fraction r\ m . We set 
the coupling strength a = 0.4£p- 

for an adiabatic sweep, but saturates at a value, that depends on initial peak phase-space 
density. Their experimental data is simulated by a Monte Carlo-type simulation 13 based on a 
model assuming that pairing formation occurs when two atoms are close enough in the phase 
space. 

On the other hand, Williams et a/. 14 discussed the molecular conversion efficiency assuming 
the molecular formation rate vanishes when the detuning 5 is negative. Their result is derived 
by using the coupled atom-molecule Boltzmann equations. They calculated the conversion 
efficiency from the molecular fraction at the detuning 5 = as a function of the initial peak 
phase space density, and found good agreement with the experimental data although they 
treated an ideal gas mixture of Fermi atoms and Bose molecules. We note that the initial 
temperature and the initial phase space density have a one-to-one correspondence. 

In the experiment for 40 K, the conversion fraction does not depend on the inverse sweep 
rate within a range from 640 to 2900fis/G. In this range, therefore the conversion process can 
be considered as adiabatic. 13 According to the kinetic theory of ref. 18, in this regime the 
chemical equilibration between atoms and molecules also brings about the thermal equilibra- 
tion. Assuming this adiabatic limit and that the principle controlling the maximum molecular 
conversion efficiency proposed by Williams et al. also applies when including the resonant 
interaction, we calculate the molecular conversion efficiency using rj m {5 = 0) as a function of 
initial temperature, as shown in Fig. 5 (A), where rj m is the molecular fraction 2N m /Nt t- We 
note that the Feshbach resonance of 40 K is known to be broad. 13 Although our mean-field 
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theory is only applicable to the narrow resonance case, we can see the qualitative effect of 
the resonance interaction from Fig. 5 (A). The dotted line represents the result from our 
model, while the dot-dashed line represents the ideal gas mixture of Fermi atoms and Bose 
molecules. 14 The dots represent the data of Hodby et al., 13 and the solid line represents their 
simulation result. Dashed lines represent the uncertainty of the pairing parameter in their 
simulation. As noted above, the weak coupling constant we used does not directly correspond 
to the experiment, and thus the comparison between the theory and experiment is only qual- 
itative. Nevertheless we find that our model agrees well with the trend of experimental data. 
The transition temperature T c at 5 = is 0.45Tp, which corresponds to the transition tem- 
perature in terms of initial temperature Tj jC = 0.36Tp. Above Tj jC our model is identical to the 
ideal gas mixture case, since we are treating the resonant interaction in a mean-field theory. 
For reference, we show the molecular fraction r] m plotted as a function of temperature and 
the detuning in Fig. 5 (B). The fraction at the resonance position 5 = and a temperature 
adiabatically connected with the initial temperature are used for the conversion efficiency in 
Fig. 5 (A). 





Fig. 5. (Color online) (A) The molecular conversion efficiencies are plotted against the initial tem- 
perature. The dotted line represents the our model, while the dot-dashed line is the result for an 
ideal gas mixture. 14 The dots represent the experimental data of Hodby et al. 13 The solid line 
represents their simulation result, and the dashed lines give the uncertainty of the pairing param- 
eter in their simulation. (B) For reference, the molecular fraction r\ m is plotted as a function of 
temperature and the detuning. The fraction at the resonance position <5 = and a temperature 
adiabatically connected with the initial temperature are used for the conversion efficiency in Fig. 
(A). We set the coupling strength a = 0.4ep- 



The resonant interaction is more important for lower initial temperatures. At the initial 
temperature T, = 0, the molecular conversion efficiency is 100% in the ideal gas mixture 
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model. On the other hand, the experimental data at Tj — > is suppressed from 100%. The 
efficiency in our model at Tj = is also less than 100%, because the resonant interaction 
suppresses the molecular conversion. We note that the resonance position may shift from the 
ideal gas case due to the atom-molecule interaction, which may also shift the result for the 
conversion efficiency quantitatively. However, the qualitative relation between the molecule 
conversion efficiency and the initial temperature will remain unchanged. 

5. Conclusion 

In this paper, we have calculated the phase diagrams for resonantly-coupled system of 
Fermi atoms and Bose molecules. We determined the paths of constant entropy traversed 
in the phase diagram as the detuning is lowered adiabatically. The adiabatic phase diagram 
of the fermionic condensate fraction composed of the condensed molecules and the Cooper 
pairs was plotted against the detuning 5 and the initial temperature of a pure atomic gas. 
The adiabatic phase diagram allowed us to compare the theory with experiments, since in the 
experients, 1,3 the condensate fraction have been measured against the initial temperature and 
the magnetic field. We obtained the transition temperature in terms of the initial temperature 
Tj jC ~ 0.21Tp in the BEC regime. Finally, the molecular conversion efficiency was plotted as 
a function of the initial temperature. We found that the resonant interaction suppresses the 
complete conversion. 

The present work is an extension of the recent works by Williams et a/. 8 ' 14 to include 
the resonant interaction using a mean-field theory. Although our model assumes a narrow 
resonance, the mean-field calculations help qualitative understanding of the recent experiments 
on the BCS-BEC crossover. In particular, the suppression of Tj jC in the BEC limit and the 
behavior of the molecular conversion efficiency will be qualitatively unchanged even in the 
broad resonance case. In order to compare the theory with experiments in a quantitative 
fashion, we must include the Cooper pairs outside the condensate, which are neglected in our 
mean-field calculation. We must also express the resonance position in terms of the magnetic 
field rather than the detuning. 
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